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ABSTRACT 

Quantum corrections to the classical field equations, induced by a scale dependent gravitatio- 
nal constant, are analyzed in the case of the static isotropic metric. The requirement of general 
covariance for the resulting non-local effective field equations puts severe restrictions on the nature 
of the solutions that can be obtained. In general the existence of vacuum solutions to the effective 
field equations restricts the value of the gravitational scaling exponent z/^^ to be a positive integer 
greater than one. We give further arguments suggesting that in fact only for z>^^ = 3 consistent 
solutions seem to exist in four dimensions. 
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Over the last few years evidence has been increasing to suggest that quantum gravitation, even 
though plagued by uncontrollable divergences in standard weak coupling perturbation theory [1], 
might actually make sense, and lead to testable predictions at the non-perturbativc level. These 
new results in general arise from the non-trivial scaling properties of the gravitational coupling 
constants in the vicinity of a non-trivial ultraviolet fixed point in four dimensions. As is often the 
case in physics, the best arguments do not come from often incomplete and partial results in a single 
model, but more appropriately from the level of consistency that various, often quite unrelated, 
field theoretic approaches provide. 

The main aspect we wish to investigate in this paper is the nature of the specific predictions 
about the running of Newton's constant G, as they apply to the standard static isotropic metric. 
Our starting point will be the solution of the non-rclativistic Poisson equation, which for a localized 
point source can be investigated for various values of the gravitational scaling exponent u. But a 
more appropriate setting will be a relativistic, generally covariant framework, wherein the effects of 
the leading quantum correction can be studied systematically, and for which we will show that the 
existence of vacuum solutions severely restricts the possible values for the exponent v. Specifically, 
we will show that no consistent solution to the effective non-local field equations can be found 
unless is an integer greater than one. To check the overall consistency of the results, a different 
approach to the solution of the covariant effective field equations for the static isotropic metric will 
be pursued, in terms of an effective vacuum density and pressure. In this case one finds that unless 
the exponent v is equal to 1/3, a consistent solution cannot be obtained. 

The starting point for our discussion is the form of the running gravitational coupling in the 
vicinity of the ultraviolet fixed point at Gc, as obtained from the lattice theory of gravity, and given 



with m = 1/^, ao > and ~ 1/3. Usually the quantity Gc in the above expression is identified 
with the laboratory scale value, ~ \fGphys ~ 1-6 x 10~^^ cm, the reason being that the scale 
^ can be very large, roughly of the same order as the scaled cosmological constant A. Quantum 
corrections on the r.h.s are therefore quite small as long as A;^ S> m^, which in real space corresponds 
to the "short distance" regime r <^^.^ For more details the reader is referred to the recent papers 
[3-5], and further references therein. 

^The result of Eq. (1) is in fact quite similar what one finds for gravity in 2 + e dimensions [6, 7, 8], 
if one allows for a different value of exponent v as one transitions from two to four dimensions, G{k^) ~ 

. See also the recent results discussed in [9, 10]. 
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For /c^ ^ the quantum correction proportional to ao diverges, and the infrared divergence 
needs to be regulated. A natural infrared regulator exists in the form of m = 1/^, and therefore a 
properly infrared regulated version of the previous expression is 

i2 
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with m = 1/^ the (tiny) infrared cutoff. Thus the gravitational coupling approaches the finite value 
Goo = (1 + ao + • • •) Gc, independent of m = 1/^, at very large distances r ^ ^. In this work we will 
be concerned with the static limit, where the non-relativistic Newtonian potential can be defined 
as 
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The static potential 0(r) can be obtained from Eq. (2) directly by Fourier transform, or equivalently 
from the solution of Poisson's equation with a point source at the origin. In the limit of weak fields, 
the relativistic field equations give for the (j) field (with goo{x) c± — (1 + 2(^{x) ) 

A(j){x) = 4TrGp{x) (4) 

and for a point source at the origin the first term on the r.h.s is just AttMG S^'^^x). The solution 
for (p(r), obtained by Fourier transforming back to real space Eq. (3), gives in the large r limit 

MG 
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with Q = l/{u22^ ^(j!/))- "^^^ part in G(k ) proportional to can equivalently be represented as 
a source term in Poisson's equation, the latter determined from the inverse Fourier transform 
of the correction term in Eq. (2), 
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One finds 



Pm{r) = ^c^aoMm^imr) -fs:i(3_ i)( 
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with Cj. = 22(5-i;)/0Fr(^). Note that the vacuum polarization density Pm(^) has the normal- 
ization property 

roo 

4:77 / r^drpmir) = cqM (8) 
Jo 

and that pm{r) diverges at small r for u > 1/3. In the small r limit and for general u > ^, one 
then finds from Poisson's equation, using the expansion of the modified Bessel function Kn{x) for 
small arguments, 

(f){r) ~ -^^^ + aoMGcsm^ r^~^ + ... (9) 



with c, = i/|sec(^)|/r(i). 

Solutions to Poisson's equation with a running G provide some useful insights into the struc- 
ture of quantum corrections, but a complete analysis requires a study of the full rclativistic field 
equations, which will be discussed next. A set of effective field equations incorporating the running 
of G is obtained from the replacement [2] 



G G{u) = G 



1 + ao - + 



^G{l + A{U)) (10) 



with the d'Alembertian □ expressing the running of G as in either Eqs. (1) or (2). The non-local 
effective field equations then read 

R^,^. - \gt,uR + \g„v = ^T^G (l + ^(n)) T^. (11) 

with ^(n) given by Eq. (10), and A ~ 1/C^- The use of the d'Alembertian □ to describe the running 
of couplings in gauge theories and quantum gravity was discussed in some detail, for example, in 
[11]. The corresponding trace equation is 

R- A\ = -SttG {1 + A{a)) T (12) 

Being manifestly covariant, these expressions at least satisfy some of the requirements for a set of 
consistent field equations incorporating the running of G. The d'Alembertian □ operator is defined 
here through the appropriate combination of covariant derivatives 

□ = g'^'V^V, (13) 

and its explicit form depends on the specific tensor nature of the object it is acting on. In general the 
operator A{D) has to be defined by a suitable analytic continuation from positive integer powers, 
which is usually done by computing for positive integer n, and then analytically continuing to 
n ^ —l/2u. Let us set for now the cosmological constant A = 0, since its contribution can always 
be added at a later stage. As long as one is interested in static isotropic solutions, one can take for 
the metric the most general form 

ds'^ = -B{r)dt^ + A{r)dr'^ + r"^ {dO'^ + sin^^V) (14) 

For the energy momentum tensor we will take the perfect fluid form 

T^v = d\&g[B{r) p{r), A{r)p{r), r"^ p{r), r"^ SIT? 9 p{r)\ (15) 
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and a point source as the origin is simply represented as 

r^,(r) = diag[S(r)/9(r), 0, 0, 0] (16) 

with the source proportional to a 3 — d delta function. 
Consider first the trace equation 

R = -87rG(n)T = -SttG (1 + A{a)) T = +87rG (1 + A{a)) p (17) 

where we have used the fact that the point source at the origin is described just by the density 
term. One then computes the repeated action of the invariant d'Alembertian on T, 

o (-8.Gr, ^ o(.8.G„ ^ 1^ - . 1^ . 5^ (18, 

In view of the rapidly escalating complexity of the problem, it seems sensible to expand around the 
Schwarzschild solution, and set 

, 1 , 2MG a(r) , , 2MG e(r) , , 

AirY^ = 1 + — ^ Bir) = 1 + (19) 

r r r r 

where the correction to the standard solution are parametrized here by the two functions cr(r) and 
0{r), both assumed to be "small", i.e. proportional to oq as in Eq. (10), with ao considered a small 
parameter. To simplify the problem even further, we will assume that for 2MG <^ r <^ ^ (the 
"physical" regime) one can set 

a{r) = -aoMGcar'^ e{r) = -aoMGcer^ (20) 

This assumption is in part justified by the form of the non-relativistic correction of Eqs. (9). Then 
for a = /3 (the equations seem impossible to satisfy if a and (3 are different) one obtains for the 
scalar curvature 

R = + a {2ca + {a - 1) ce) ao M G r""^ + +0(ag) (21) 

A first result can be obtained in the following way. Since in the ordinary Einstein case one has for 
a perfect fluid R = —SirGT = +87rG(p — 3p), and since Pm{^) ~ r^~^ from Eq. (7) in the same 
regime, one concludes that a solution is given by 

a = ^ (22) 

which is also consistent with the Poisson equation result of Eq. (9). 

The next step up would be the consideration of the action of □ on the point source, as it appears 
in the full effective field equations of Eq. (11), with again T^i, described by Eq. (16). One perhaps 



surprising fact is the generation of an effective pressure term by the action of □, suggesting that 
both terms should arise in the correct description of vacuum polarization effects, 

pB'"^ 2Bp' BA'p' B' p' Bp" 

(□T;.),, = (23) 

and (□ T^^)gg = iOT^v),^^p = 0. A similar effect, namely the generation of an effective vacuum 
pressure term in the field equations by the action of □, was seen already in the case of the Robertson- 
Walker [2]. 

To check the overall consistency of the approach, consider next the set of effective field equations 
that are obtained when the operator (1 + A{B)) appearing in Eqs. (11) and (12) is moved over 
to the gravitational side. Since the r.h.s of the field equations then vanishes for r 7^ 0, one has 
apparently reduced the problem to one of finding vacuum solutions of a modified, non-local field 
equation. Let us first look at the relatively simple trace equation. If we denote by SR the lowest 
order variation (that is, of order ao) in the scalar curvature over the ordinary vacuum solution 
R = 0, then one has 

On a generic scalar function F{r) one has the following action of the covariant d'Alembertian □ : 

A'F' B'F' 2F' F" , , 

°nO = -^ + ^ + ^ + ^ (25) 

Assuming a power law correction, as in Eq. (20), with a = /3, as in Eq. (21), one then finds 
UR ^ aoMGr"-^ (2ca + C0(a-l)) a(a-2)(a-3) 

n^R aoMGr"-'^ (2C(, + ce(a- 1)) a(a-2)(a-3)(a-4)(a-5) (26) 

and so on, and for general ^ — ^ +^ 

U-R ^ aoMG(2c, + c,(a-l)) ^^^g,ii^r«-^-=^ (27) 

V\l — a) 

Therefore the only possible power solution for r ^ MG is a = 0, 2. . . ^ + 1, with c„ and cq 
unconstrained to this order. 

Next we examine the full effective field equations (as opposed to just their trace part) as in 
Eq. (11) with A = 0, If one denotes by 50^^ = 6 (yR^v — 5 Qixv R^ the lowest order variation (that 
is, of order ao) in the Einstein tensor over the ordinary vacuum solution Gfxv = 0, then one has 

=0 (28) 



again for r 7^ 0. Here the covariant d'Alembertian operator □ acts on a second rank tensor and 
would thus seem to require the calculation of as many as 1920 terms, of which many fortunately 
vanish by symmetry. In the static isotropic case the components of the Einstein tensor are given 

by 

^ _ A'B B B 
G = -- + — + - 

1^2 y> 1^2 

B'\^ A'B'r^ B"r^ A'r B'r 

Gee = -TT7^- . .9^ +77T7T-7r7^ + 



4^425 2AB 2A^ 2AB 



G^^ = sin^eGee (29) 



After acting with □ on this expression one finds a rather complicated result. Here we will list only 

6BA'^ 2BA'^ AB'A'^ 2BA' 6BA"A' B"A' 6B 
+ ^r-n 71 ^ 7T— + — 7^ + 



rA^ r'^A* rA'^ r^A^ rA'^ rA^ r^A 
6B _ AB^ 4B^ _ BA^ 2B'A" B" _ BA^^) 

A A^ T T^A T^A^ tA^ 

If one again assumes that the corrections are given by a power, as in Eq. (20), with a = (3, then 
one has to lowest order 

Gtt = aoMGc^ar'^-^ 

Grr = -aQMG{c^+C0{a-l))r''-^ 

Gee = -^aoMG(c^ + ce(a-l))(Q-l)r"-i (31) 

with the ipip component again proportional to the 99 component. Applying □ on the above Einstein 
tensor one then gets 

(□G)^j = aoGMc^a(a-2)(a-3)r"-^ 

{nG)„ = -a^GM {ca + ce{a-l)) a{a-?>)r''-^ 

{UG)ee = -^aoGM (c^ + ce(a-l)) a(a-3)2r"-3 (32) 



(with the (fif component proportional to the 99 component), and so on. One then has for general 

r(2 + i 



(a - 1) r(-a) 
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(□"G),, ^ -aoGM ic^ + ce{a-l)) \ ^ ' r 

(a - 1) \Oi--) r(-a 



1 fa-l-i) r(2 + i-a) , 

(□"G),, ^ a,GM{c, + ce{a-l))^ \ r""^" (33) 

^ (a - 1) (^a- - j r(-a) 

Inspection of the above results reveals a common factor l/r(— a), which would allow only integer 
powers a = 0, 1, 2 . . ., but the additional factor of 1/(q! — 1) excludes a = 1 from being a solution. 
Even for a close to \/v (as expected on the basis of the non-relativistic expression of Eq. (9), as 
well as from Eq. (22)) ~ 1/a — e only integer values a = 2,3,4. . . are allowed. In general the 
problem of finding a complete general solution to the effective field equations by this method lies 
in the difficulty of computing arbitrarily high powers of □ on general functions such as C7(r) and 
0(r), which eventually involve a large number of derivatives. Assuming for these functions a power 
law dependence on r simplifies the problem considerably, but also restricts the kind of solutions 
that one is likely to find. More specifically, if the solution involves (say for small r, but still with 
r » 2MG) a term of the type In mr, as in Eqs. (9), (48) and (51) for v 1/3, then this 
method will have to be dealt with very carefully. This is presumably the reason why in some of 
the F-function coefficients encountered here one finds a power solution (in fact a = 3) for u close 
to a third, but one gets indeterminate expression if one sets exactly a = 1/z^ = 3. 

The earlier discussion of the non-relativistic case suggests that the quantum correction due 
to the running of G can be approximately described by Poisson's equation, with a source term 
related to a vacuum energy density /5m(r), distributed around the static source of strength M in 
accordance with the result of Eqs. (7) and Eq. (8). These expressions, in turn, were obtained by 
Fourier transforming back to real space the original result for G(fc^) of Eq. (2). Furthermore, in 
the preceding discussion of the relativistic case it was found (as in [2] for the Robertson- Walker 
metric case) that a manifestly covariant implementation of the running of G, via the G(n) given in 
Eq. (10), will induce a non- vanishing effective pressure term in the field equations. This result can 
be seen clearly, in the case of the static isotropic metric, for example from the result of Eq. (23). 
We will therefore now consider a relativistic perfect fluid, with energy-momentum tensor, which in 
the static isotropic case reduces to Eq. (15). The tt, rr and 06 components of the field equations 
then read 

, „ A'B B B 
-\B + —-- — + — = SttGBp 



A B' 1 
B'\^ A'B'r^ B"r^ A'r B'r 



with the iptp component equal to sin^^ times the 06 component. Energy conservation V^T^j/ = 



implies 

(P + P)^+P' = (35) 

and forces a definite relationship between B{r), p{r) and p{r). The three field equations and the 
energy conservation equation are, as usual, not independent, because of the Bianchi identity. 

It seems reasonable to attempt to solve the above equations (usually considered in the context 
of relativistic stellar structure) with the density given by the Pm{f) of Eqs. (7). This of course raises 
the question of how the relativistic pressure p{r) should be chosen, an issue that the non-relativistic 
calculation did not have to address. We will argue below that covariant energy conservation com- 
pletely determines the pressure in the static case, leading to consistent equations and solutions 
(note that in particular it would not be consistent to take p{r) = 0). Since the function B[r) drops 
out of the tt field equation, the latter can be integrated immediately, giving 

r 3 r Jo 

It also seems natural here to identify ci = —2MG, which of course corresponds to the correct 

solution for qq = {p = p = 0). Next, the rr field equation can be solved for B{r), 

W ^ I r ^ 1 + A(y) {\y^-8nGy^p{y)-l) \ 

B(r) = exp C2 - / dy ^ '- ) (37) 

\ Jro y J 

with the constant C2 again determined by the requirement that the above expression for B{r) reduce 
to the standard Schwarzschild solution for oq = (p = p = 0), giving C2 = ln(l — 2MG/ro — Xtq/S). 
The last task left therefore is the determination of the pressure p{r). Using the rr field equation, 
B'{r)/B{r) can be expressed in term of A{r) in the energy conservation equation. Inserting then 
the explicit expression for A{r), from Eq. (36), one obtains 



U^Gr^pir) + 2MG - \\r^ + SttG T dxx'^p{:x)] {p{r) + pir)) 

p'{r) + ^ ^ °- ^ = (38) 

2 r (r - 2MG - | - SttG dx x'^ p(x)) 

which is usually referred to as the equation of hydrostatic equilibrium. From now on we will focus 

only the case A = 0. The last equation, a non-linear differential equation for p{r\ can be solved to 

give the desired solution p(r), which then, by equation Eq. (37), determines the remaining function 

B{r). In our case though it will be sufficient to solve the above equation for small oq, where oq (see 

Eqs. (2) and (7)) is the dimensionless parameter which, when set to zero, makes the solution revert 

back to the classical one. It will also be convenient to pull out of A(r) and B(r) the Schwarzschild 

solution part, by introducing the small corrections air) and ^(r), as defined in Eq. (19), both of 

which are expected to be proportional to the parameter oq. One then has 



Again, the integration constant C2 needs to be chosen here so that the normal Schwarzschild solution 
is recovered (oi p = p = 0. To order ao the resulting equation for p{r), from Eq. (38), is 

Note that in regions where p{r) is slowly varying, p'{r) ~ 0, one has p 2± —p, i.e. the fluid 
contribution is acting like a cosmological constant term with cr(r) ~ 6{r) ~ —{p/3)r^. The last 
differential equation can then be solved for p{r) , 

= L - r _M£^] (41, 

where the constant of integration has to be chosen so that when p(r) = (no quantum correction) 
one has p{r) = as well. Because of the singularity in the integrand at r = IMG, we will take 
the lower limit in the integral to be rg = 2MG + e, with e ^ 0. To proceed further, one needs the 
explicit form for pm{r), which was given in Eqs. (7) and (8). The required integrands involve for 
general v the modified Bcsscl function Kn{x), and can be therefore a bit complicated. Here we will 
limit our investigation to the small r {mr <S 1) and large r {mr S> 1) behavior. Since m = 1/^ is 
very small, the first limit appears to be of greater physical interest. For small r the density Pmif) 
has the following behavior (see Eq. (7)), 

Pm{r) - (42) 

r— >0 

for 1/ > 1/3, with the constant 

I sec f—) I 1 

^0 = aoMmi^ (43) 

47rr(i-l) 

determined from the small x behavior of the modified Bessel function Kn{x). For u <\/3 pm{r) ~ 
const. aoMm^, independent of r. For = 1/3 the expression for pm{r) is given later in Eq. (46). 
Therefore in this limit, with ^ < v < 1, one has 

A-^{r) = 1 - - 2aQMGcsm^ r^-^ + ... (44) 

with the constant Cg = i'\ sec (^) |/r(i — 1). For i/ = 1/3 the last contribution is indistinguishable 
from a cosmological constant term — |r^, except for the fact that the coefficient here is quite 
different, being proportional to ~ aoMGm^. To determine the pressure, we suppose that it as 
well has a power dependence on r in the regime under consideration, Pm{f) = Cp Aq r'^, where Cp is 
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a numerical constant, and then substitute Pm{^) into the pressure equation Eq. (40). This gives, 
past the horizon r 3> 2MG, 

(27 - l)cpMGr^-^ - Cpjr'^ - MGr^/"'^ ~0 (45) 

giving the same power 7 = l/i/ — 3 as for p{r), Cp = —1 and surprisingly also 7 = 0, implying that in 
this regime only 1/ = 1/3 gives a consistent solution. Again, the resulting correction is quite similar 
to what one would expect from a cosmological term, with an effective Xm/-^ — STTuaoMGrn^. 
The case = 1/3 requires special treatment, and one needs to go back to the expression for 

Pm{r) for u = 1/3, 

1 -i 

Pm{r) = ^aoMm-^ Ko{mr) (46) 

For small r one then has 



Pmir) = -^Mm' (In— + 7 j + ••• (47) 

and consequently from Eq. (36), 

A-^{r) = 1 + — ^ In (mr) + ... (48) 

r 3 TT 

Prom Eq. (40) one then obtains an expression for the pressure Pmir), and one finds 

3i ^ ^ aoMnr" log f r + r Jl - - MG] ,,3 ,.3 
_ aQMmr log(mr) ^ ^ \^ ' V r j aQMmr aoMnv^cz 

Pm y) — n o / = I o I / = 

27r2 27ryi - ^ 27ryi - ^ 

'(49) 

where C3 is again an integration constant. Here we will be content with the r 3> 2MG limit of the 
above expression, which we shall write therefore as 



PmiT) = — —^Mm^ In (mr) + ... (50) 

After performing the required r integral in Eq. (39) , and evaluating the resulting expression in the 
limit r 3> 2MG, one obtains an expression for 9{r), and from it 

2MG AaoMGm^ 21/ x 

B{r) = 1 1 ^ In (mr) + ... (51) 

r 3 vr 

The expressions for A{r) and B{r) are, for r S> 2 MG, consistent with a gradual slow increase in 
G in accordance with the formula 

G ^ G{r) = G (1 + ^m^r^ In + ..) (52) 
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and therefore consistent as well with the original result of Eqs. (1) or (2), namely that the classical 
laboratory value of G is obtained for r ^ ^. In fact it is quite reassuring that the renormalization 
properties of G{r) as inferred from A{r) are the same as what one finds from B{r). For large r one 
has instead, from Eq. (7) for pm{r), 

Pm{r) - Aor^-^e-"^' (53) 

r — > oo 

with ^0 = l/V 1287r ao M m}'^^ . In the same limit, the integration constants is chosen so that 
the solution for A{r) and B{r) at large r corresponds to a mass M' = (l + ao)M (see the expression 
for the integrated density in Eq. (8)), or equivalently 

air) ~ eir) ~ -2aoMG (54) 

r ^ oo 

One then recovers a result similar to the non-relativistic expression of Eq. (5), with G{r) approach- 
ing the constant value G^o = (1 + clq) G, up to an exponentially small correction in mr at large 
r. 

In conclusion, it appears that a solution to the relativistic static isotropic problem of the running 
gravitational constant can be found, provided that the exponent u in either Eq. (2) or Eq. (11) 
is close to one third. This last result seems to be linked with the fact that the running coupling 
term acts in some way like a local cosmological constant term, for which the r dependence of the 

vacuum solution for small r is fixed by the nature of the Schwarzschild solution with a cosmological 
constant term. Furthermore, in d > 4 dimensions the Schwarzschild solution to Einstein gravity 
with a cosmological term is given by [12] 

{d-2)iT— (d-2)(d-l) 
which would suggest, in analogy with the results for d = 4 given previously, that in d > 4 dimensions 
only 1/ = l/(d — 1) is possible, if the correction again behaves locally like a cosmological constant 
term. This last result would also be in agreement with the exact value v = found at d = oo [5]. 

To summarize, the starting point for our discussion of the renormalization group running of G is 
Eq. (1), valid at short distances k ^ m, oi its improved infrared regulated version of Eq. (2). While 
a solution to the non-rclativistic Poisson equation can be given for various values of the exponent v, 
the scale dependence of G can also be consistently embedded in a relativistic covariant framework 
using the d'Alembertian □ operator. This then leads to a set of nonlocal effective field equations, 
whose consequences can be worked out for the static isotropic metric, at least in a regime where 
2MG ^ r ^ ^, and under the assumption of a power law correction. We have found that the 
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structure of the leading quantum correction is such that it severely restricts the possible values for 
the exponent u, in the sense that no consistent solution to the effective non-local field equations, 
incorporating the running of G, can be found unless is an integer. A somewhat different 
approach to the solution of the static isotropic metric was pursued in terms of an effective vacuum 
density of Eq. (7) , and a vacuum pressure chosen so as to satisfy a covariant energy conservation for 
the vacuum polarization contribution. The main result there is the derivation, from the relativistic 
field equations, of an expression for the metric coefficients A{r) and S(r), given for 2MG <^r 
in Eqs. (48), (51) and (52). Prom the nature of the solution for A{r) and B{r) one finds again that 
unless the exponent v is close to 1/3, a consistent solution to the field equations cannot be found. 
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